Stark-modulation spectroscopy of the B(1)[ 3 II] state of PbO 
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We report detailed spectroscopic measurements of the X(0) (v = 0) — * B(l) [ 3 Il] (v = 5) 

transition in PbO. Using a Stark-modulated laser absorption technique, we have measured the 
hyperfine constant of 207 PbO in the B(l) state, as well as the B(l)(v — 5) rotational constant, 
X — B isotope shifts, etc. The hyperfine constant of the B(l) state is of interest as a benchmark 
for calculations of PbO electronic structure, related to experiments to search for the electric dipole 
moment of the electron. 

PACS numbers: 33.15.-e,33.20.Kf,14.60.Cd 



Spectroscopy of the PbO molecule has become of in- 
terest since it is a good candidate for use in searches for 
a permanent electric dipole moment (EDM) of the elec- 
tron, d e . EDMs are interesting because a non-zero value 
for an EDM of any fundamental particle would violate 
parity and time-reversal symmetries. Standard model 
predictions for EDMs are well below any proposed ex- 
perimental sensitivities. However, most extensions of the 
standard model predict dramatically enhanced EDMs. 

A promising approach towards improving the limits on 
d e involves heavy polar diatomic molecules in a configu- 
ration with an unpaired electron spin. In such systems, 
the effective electric field Wd seen by the unpaired elec- 
tron can be many orders of magnitude larger than exter- 
nal fields attainable in the laboratory 0, Q ■ An electron 
EDM can be detected by spin polarizing an electron along 
this internal field and searching for the characteristic lin- 
ear Stark shift AE — —d e Wd- Interpreting experimental 
limits on AE measured in a molecule in terms of d e re- 
quires a value for Wd- This can be obtained using semi- 
empirical wavefunctions for the state of interest pj, or 
from ab initio calculations 0- Spectroscopic properties 
sensitive to the electron spin density at the nucleus such 
as hyperfine structure (hfs) can constrain the parameters 
in semi-empirical evaluations of Wd, or test the predic- 
tions of ab initio calculations. 

An experiment is underway in our lab using the 
o(l) [ 3 S + ] state of PbO, that may improve the current 
limit I d e | < 1.6 x 10 -27 e cm 0] by several orders of 
magnitude Such an improvement would probe large 
regions of the parameter space of many standard model 
extensions. Stimulated by this work, a measurement of 
the hyperfine constant of the a(l) state was recently re- 
ported @. The hfs of the B(l) [ 3 IJ] state of 207 PbO 
is also of interest, since it can be used to further check 
and/or refine the electronic structure calculations used to 
calculate the value of Wd in the a(l) state. In addition, 
knowledge of the B(l) hfs can be used for estimating Wd 
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for the B(\) state, which also may be a viable candi- 
date for an EDM search ||| . This motivates this effort to 
determine the hfs of the B(l) state of PbO. 

The measurements were made by observing the ab- 
sorption of laser light by PbO vapor in a cell as the laser 
wavelength was tuned. The cell contains PbO of natural 
isotopic abundance; it is constructed from a 5.7 cm cube 
of alumina, with 5.1 cm diameter through holes bored 
perpendicular to, and centered on, each face. The top 
and bottom holes of the resulting cubical frame struc- 
ture are capped with alumina plates. Thin YAG win- 
dows are bonded on the remaining four holes using gold 
foil as an intermediate layer. The cell is heated to 700°C 
by the radiation from resistively heated thin foils of tan- 
talum supported on a quartz framework surrounding the 
cell. This results in a PbO number density n«4x 10 13 
cm~ 3 in the cell. The laser used in the measurements 
was an external cavity diode system based on a standard 
design Q. Our system used a Nichia NLHV500C violet 
diode laser and 3600 grooves/mm diffraction grating in 
a Littrow configuration. The laser provided <1 mW of 
power around 406.5 nm which passed through the YAG 
windows into the cell. 

To enable sensitive lock-in detection of the absorption, 
we employed a Stark modulation technique As ex- 
plained below, this method enhances signals from low J 
levels where the hfs is largest, which are otherwise diffi- 
cult to isolate within the spectrally congested region near 
the bandhead. Here we briefly discuss the basic features 
of the Stark- modulated signals, as relevant to our exper- 
iment. 

The Stark effect in the X(0) — > B(l) transition is dom- 
inated by shifts due to mixing between the closely-spaced 
O-doublet levels in the B(l) state. The doublet mem- 
bers, e and /, share a common value of total (electronic 
+ rotational) angular momentum J but are of opposite 
parity. The spacing between the doublet levels is given 
by Acj(J) = qJ(J + 1), where q <C B v=5 , the rotational 
constant for the B{v — 5) state. [By convention, q is pos- 
itive if the e level, with parity P = {— 1) J , is at higher 
energy.] In the absence of an electric field E, selection 
rules for the electric-dipole X — > B transition ensure that 
the laser light is absorbed only in the transition to one 
of the doublet levels. In the presence of E, the e and / 
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FIG. 1: Stark modulation of the excited state energy levels, 
leading to a modulation of the transmitted light intensity. 
Plot (A) shows the Ql transition and the absorption expected 
at E=0. In (B), with E large, the Stark effect mixes and shifts 
the doublet levels. The difference in absorption is recorded in 
the scans. 



levels mix and are shifted by an amount Av given by 
Av = ± 



V /(A n (J)/2) 2 + (^)2 - A (J)/2 



(1) 



Here the upper (lower) sign corresponds to the upper 
(lower) level in the limit E — > 0, and 



a = Hb(i)M f 



F{F + l) + J{J+l) -1(1+1) 
2F(F + 1)J(J + 1) 



(2) 



where /j,b(i) is the dipole moment of the B state, and 
F = J + / is the total angular momentum. In the limit 
of large E and small J, Av oc E and the doublet levels 
are completely mixed. In this limit, the laser is absorbed 
with equal strength by both eigenstates of the mixed dou- 
blet, and the repulsion of the doublet levels results in a 
characteristic broadening of the absorption spectral line 
(see Fig. Q). For a laser tuned to one of the unperturbed 
doublet levels, by modulating E, maxima in absorption 
will occur twice per period at the minima of E. The f2- 
doublet mixing leads to small shifts in the line centers, as 
does the slight residual Stark mixing with nearby rota- 
tional levels of opposite parity. Both of these small shifts 
(< 50 MHz) are accounted for in the analysis. 

In our measurements, an electric field E was produced 
using two parallel 3.8 cm diameter gold foil electrodes 
spaced 0.6 cm apart in the cell. The laser light was lin- 
early polarized parallel to E and traversed the vapor cell 



between the electrodes. A reference signal at u> = 2ir x 99 
kHz was amplified to produce a voltage V — Vq cos(ujt), 
where Vb ~ 800 V, across the electrodes. The transmit- 
ted light was detected with a photodiode and amplified. 
Part of the laser beam was diverted for a wavelength 
measurement (using a Burleigh WA-1500 Wavemeter), 
and the laser power reflected from the cell was monitored 
with a second "reference" photodiode. 

The transmitted signals were sampled at 1 MHz with 
a 12 bit ADC and recorded for 4 s. The amplitude of the 
signal at the 2nd harmonic of the reference was extracted 
from the digitized data. Then the laser frequency was ad- 
vanced roughly 150 MHz by simultaneous adjustment of 
the grating angle and diode current, and a new absorp- 
tion measurement was taken. This cycle was repeated 
until a mode- hop occurred, yielding scans covering 15-30 
GHz. 

The scans showed absorption features correspond- 
ing to 208 PbO, 207 PbO, and 206 PbO molecules making 
P-, Q-, and i?-branch transitions from the electronic 
ground state X(0) [ X S+] (v" = 0, J") to the excited state 
5(1) [ 3 n] (v' = 5, J'), with S/N > 1 for J" = — 6. An 
example is shown in Fig. [21 Before analyzing the data, 
the absorption amplitudes were corrected for the varia- 
tion in laser power with frequency during a scan. We 
expect residual variations in amplitude over the scan of 
~ 5%, based on observed fluctuations in the ratio be- 
tween transmitted and reflected power (presumably due 
to interference effects in the cell windows). The statisti- 
cal uncertainty in each measurement was estimated from 
the fluctuations of the second harmonic out of phase with 
the signal, and was typically a few times greater than 
the shot noise in the photodiode signal. The largest sig- 
nal was observed from the Ql 208 PbO line, with a peak 
S/N >75, corresponding to » 4 ppm absorption. These 
scans were fit to extract the hfs constant and coefficients 
in the Dunham expansion, as described below. 

Hyperfine structure in the B(l) state of 207 PbO arises 
from the interaction Hm s = J c • A • I = A\\ J c zlz + 
A± (Jcxlx + JcyW), where A is the hyperfine tensor, 
I = 1/2 is the nuclear spin of 207 Pb, J e is the electronic 
angular momentum, and the coordinate axes are defined 
in the body-fixed frame, with Z along the internuclear 
axis n. For levels with low J values in an £1 = 1 elec- 
tronic state, the hfs is dominated by the term in -f/hfs 
proportional to Au. The terms containing Jx,y induce 
off-diagonal mixing with states with f2 ^ 1, which are 
far in energy for this Hund's case (c) molecular state. 
We can thus write the energy shift due to hfs, correct to 
second order in H^s ~ A\\J e zIz> as 
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The matrix elements are determined from: 
(FIJ'M F fl\J eZ Iz\FIJM F n) 

\F+/+J (F I J' 



i j i 



y/l(I+l)(2I + l) 



x(-i)^-V(2J' + i)(2J + i)( J fi I |Hfi. (4) 



Eqn. |3 then simplifies to: 



TABLE I: Dunham coefficients V. 



B.20S 



in cm . The value 



given in the place for Yoo is actually T e + Yim 



i \ J 





1 





22282.4297(60) 


0.26465(1) 


1 


498.497(13) 


-0.002561(12) 


2 


-2.331(63) 


-3.4(4) xl0~ 5 


3 


0.03638(80) 





A<r ±1/2 = ± 



2(F + l/2) 

(F-l/2)(F + 3/2) 



&B V = 5 



(F + 1/2)3 



(5) 



Initial assignments of each line were possible using es- 
timates of the B(l) state energies by extrapolation from 
the data of Ref. 9] , and taking into account the known 
natural abundances of each Pb isotope. The line splitting 
due to hfs in 207 PbO is clearly visible in two transitions, 
corresponding to the Ql and Q2 lines, in FigEl T ne un- 
perturbed line positions of these hfs-split lines could be 
located initially by extrapolating from 207 PbO lines with 
higher values of J, which exhibit much smaller hfs split- 
ting. The sign of Au is found to be positive by visual 
inspection of Fig. [2| using the following observations: 
the largest peak, which is closest to the unperturbed line 
center, must arise from the F = J +1/2 state (see Eqn|S}; 
and in addition, if An > (< 0), then this line should 
lie higher (lower) in energy than the unperturbed line 
center. 

Once the lines were assigned, each scan was analyzed 
with a multiparameter fit in which the modulated ab- 
sorption lines were approximated as as sum over second 
derivatives of a gaussian, with a separate gaussian for 
each value of mj in the excited state. With the iso- 
topic abundances fixed, several fit parameters were used 
to characterize the height and the width of the various 
transitions; these included the ratio olP:Q and R.Q tran- 
sition strengths, the value of Hb(i), the Doppler width, 
and a measure of the electric field inhomogeneity. To ex- 
tract line positions, the energy E Z ' A (v,J) of electronic 
state Z (Z = X, B) with vibrational and rotational quan- 
tum numbers v and J, respectively, for the A PbO isotopic 
species was written as 



E Z ' A (v, J) =Gf A + B Z ' A J (J + 1) - D Z ' A J 2 (J + 1) 



E x - A {v = 0, J = 0). 



(6) 



The energies of the X state sublevels were taken from the 
precise data of Ref. |9| . The B state values were deter- 
mined by fitting, with separate values of the rotational 
constants used for the Sl-doublet levels with e and / char- 
acter. For 207 PbO, the hfs constant A\\ was included as 
an additional fit parameter. 

As described above, determination of An was the pri- 
mary goal of our work. Fitting the data set yielded 
An = 5.01(7) GHz, where the dominant uncertainty is 



from the spread in the fit results from different scans, 
with smaller contributions from variations in the fit re- 
sults from changes in the laser power correction and 
other input parameters. This value is in good agree- 
ment with the ab initio prediction in 0. We also find 
G*f4° 8 = 24600.085(8) cm" 1 , (where the uncertainty 
is from the wavemeter absolute calibration specification, 
which we did not verify independently, and from the ex- 
traction of the line centers); B^lf 8 = 0.24939(20) cm -1 ; 
and Hb(i) = 4(1) D. Our data is insensitive to the small 
parameters D^® 8 and q. 

The shifts between the values of G^l 5 for the different 
isotopes arise from changes in both the rovibrational and 
electronic energies. In order to separate these effects, it is 
necessary to evaluate the absolute ro- vibrational energies 
of the B(l)(v = 5) and the X(0)(v = 0) states. We 
accomplish this using a Dunham expansion : 

E z > A (v, J) = T Z ' A + Yf k A (v + l/2) j (J (J + l)f . 
j'.fc 

(7) 

Here the isotopic changes in rovibrational energy are 
taken into account by writing = Y^ z k ' 20S p A ^ '/ 2+fc ^ 

where p& = //(A)//i(208) and /i(A) is the molecular 
reduced mass of A PbO. The coefficients in the Dun- 
ham expansion are simply related to the previously de- 
scribed parameters. For example, — T e B ' A + 

Ej yf, A (5 + 1/2) J - Ej Y? Q A (0 + l/2)\ 

To evaluate the rovibronic part of the isotopic shifts, 
we take the values of Y^' 208 from Ref.0- The Dunham 
coefficients Y^' 208 given in Table [I] were evaluated using 
the data of Ref. 9] on the three lowest vibrational levels 
(v = 0, 1,2) of the B(l) state, in combination with our 

B 208 B 208 a 

determination of G v ^ 5 and B ' ' 5 . Although inclusion 
of our data shifts the values of individual Dunham coeffi- 
cients noticeably (as compared to their values using only 
the data from Ref. |9j), the resulting change in the ro- 
vibronic part of the isotope shifts is small: < 100 MHz, 
comparable to the uncertainty in our determination of 
individual line positions. Due to our limited range in J, 
we have negligible sensitivity to isotopic differences in the 
rotational constants. 

With this determination of the rovibronic part of the 
isotope shift, it is possible to extract the isotopic shifts 
in the electronic energies. Defining Av(A — A') = 
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Laser Frequency in GHz - 737000 GHz 

FIG. 2: Segment of a typical Stark modulation scan. The points are data and the line is a fit. The vertical bars mark the line 
centers in the absence of hfs; the hyperfine splitting of the 207 PbO lines is clearly visible. 



( t b,a_ t x,A) _ f T B,A' _ Tf> A '), we find Aj/(208 - 

207) = -270(50) MHz and A^(208 - 206) = -720(60) 
MHz. In a second approach, line centers from two differ- 
ent isotopes were extracted with no constraints on their 
relative positions or amplitudes. Subtracting off the rovi- 
bronic part of the separation expected between the lines 
from Eqn[7|gave consistent results for the electronic part 
of the isotope shifts. To account for these electronic 
shifts, we write 

T^ = T e z ' p ^l + V z (rl) + A z ^y (8) 

In this expression, T z,p is the electronic energy of state 
Z for a hypothetical point-like nucleus of infinite mass, 
and the parameters V z and A z characterize the field- 
shift and mass-shift parts, respectively, of the electronic 
isotope shift [l3l Il4|: here (r\) and Ma are the mean- 
square charge radius and the mass, respectively, of the 
A Pb nucleus, and m e is the electron mass. We use the 



known difference in mea n-sq uare charge radii between 
the different Pb isotopes [15j to separate our measured 
values for the electronic shift into field- and mass-shift 
contributions. We find V B - V x = 1.0(7) x 10~ 5 fm~ 2 , 
and A B - A x = 8.5(3.1) x 10" 6 . The normal part of 
the mass shift (due to the change in the electron reduced 
mass) is negligibly small; hence the observed shift is a 
specific mass shift, due to electron correlations. This shift 
is unusually large, although similar shifts are known to 
occur in complex atoms |14| . 

In summary, we have measured some detailed spectro- 
scopic parameters of the X(0)(v = 0) — > B(l)(v = 5) 
transition of PbO. The measured value of the B(l) state 
hfs provides an important data point for checking calcu- 
lations of the effective electric field acting on an electron 
EDM in excited states of PbO. 
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